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Abstract 

In B Xg^ there are perturbative QCD corrections at order Og not only to 
the single quark line process 6 — > 57 but also coming from a set of diagrams 
where the weak interaction vertex involves a gluon which interacts with the 
spectator quark in the B hadron. We discuss the impact of the full set of 
these diagrams. These can influence the decay rate and also the line shape 
of the photon spectrum as they favor a softer photon energy than does the 
pure spectator decay b — > 57. A subset of these diagrams generate differences 
in the decay rates for charged and neutral B's , — > Xg'j vs. 5° ^ A,7 
which can be searched for experimentally, although our calculation indicates 
that the effect will be hard to observe. Bound state effects turn out to be 
required to give rise to spectator corrections to the decay rate at order a^. 
The resulting contribution to the total rate is rather small (of about 5% ), 
although the effect is somewhat larger on portions of the photon spectrum. 
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1 Introduction 



Rare decays of heavy mesons received considerable interest in the recent years as 
possible grounds for testing the Standard Model in flavor changing neutral current 
processes which known to occur only at one loop level Exclusive decays of B- 
mesons to the non-charmed final states (which are examples of FCNC) are somewhat 
easier to observe experimentally. However theoretical calculations of exclusive rates 
include modeling of the hadronization of the final s((i)-quark and the resulting pre- 
dictions vary by the order of magnitude from one to another. The inclusive radiative 
decay B —>■ Xgj was recently observed by CLEO collaboration and is consistent 
with the SM expectations. It is now generally accepted that QCD corrections sig- 
nificantly change the decay rate for this process so it is important to bring all of 
the possible corrections under control to reduce theoretical uncertainties associated 
with QCD. 

Most work on the inclusive B X^'^ decay proceeds by calculating radiative 
corrections to the quark decay h 57 (see for example [Q). However, there are 
also several diagrams where a radiative correction involves a gluon which interacts 
with the spectator quark. They generate effects which are not accounted for in 
spectator model calculations, and the inclusive rate B Xg'y could in principle 
differ from that of 6 — S7 at order a^. If one is interested in a complete analysis of 
this decay mode, these non-spectator diagrams must be included. This is especially 
true if one is interested in the photon spectrum, since reactions which include the 
spectator quark produce a softer photon spectrum than does the free quark decay. 
The calculation presented below fills this gap and completes the calculation of the 
process of inclusive radiative decay of heavy meson at the level 0{as) 

In this paper, we present a QCD based calculation which estimates these non- 
spectator effects. There are 12 diagrams which we consider, shown in Fig 1. The 
weak vertices are calculated in QCD perturbation theory (some vertices already exist 
in the literature), and the ACCMM model |^ is used to calculate the decay rates 
and spectra. The interference of these diagrams with the leading 6 — 57 amplitude 
only occurs for small values of the momentum of the spectator quark. This portion 
is then particularly sensitive to the bound-state wavefunction, but fortunately the 
interference is relatively small. The main non-spectator effect comes from the region 
of phase space where the spectator quark recoils significantly, and hence are second 
order in a^- 

The plan of the paper is as follows. In Section 2, we organize the calculation and 
recall some of the information which in known about the gluonic vertices. Section 
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3 is devoted to most of the details about the calculation of the modifications to the 
decay rate and photon spectrum. The impact of the results is described in Section 
4. Before proceeding, let us resolve a nomenclature ambiguity. It is conventional 
to call calculations which focus on the b ^ s'j transition by the phrase "spectator 
model" and to refer to the light quark in the B hadron as the "spectator" quark. We 
will continue to use the latter phrase for the light quark in the B meson although, 
in the diagrams which we consider, it is no longer a spectator but is an active 
participant in the weak decay. We will call the 6 — > 57 diagrams the spectator 
process and conversely will label those which actively involve the spectator quark 
as non-spectator transitions. 

2 Diagrams and vertices 

There are quite a number of diagrams at order which lead to modifications to 
B — > Xs7, as shown in Fig 1. Let us divide these into five classes: 

1. QC-D-corrections to the single quark process b — > 57. 

2. "Gluonic-penguin" (bremstrahlung) diagrams which include bsg-vertex with 
photons emitted from external quark legs. 

3. "Triangle" diagrams in which both the gluon and photon are emitted from 
weak vertex due to a quark loop. 

4. Hard gluon spectator exchange diagrams with gluons emitted from b or s-lines. 

5. Diagram (10) involves the H^H^7-vertex and and is suppressed by a power of 

We will not consider this diagram further. 
The calculation proceeds by evaluating each of the diagrams, adding them to- 
gether coherently, integrating over the initial particles momentum using a momen- 
tum space wavefunction and summing over the final momenta and spins. The effec- 
tive vertices are often momentum dependent and depend on the momentum flow in 
the diagram. Since we are dealing with inclusive decays, the details of hadronization 
of the final s-quark state are not relevant for our discussion (all hadronic final states 
are summed over), so we have to model only the initial state. This substantially 
reduces the model-dependence of the final result. However we do have to carefully 
account for interference effects between the many diagrams, as different diagrams 
contribute to different portions of the final particles' phase space. For the 6 — S7 
process, the spectator quark emerges with only low momentum reflecting its mo- 
mentum in the bound state wavefunction. However in the other diagrams, the final 
quarks carry a typical momentum of a fraction of the B mass. Interference between 
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the two classes only occurs in a portion of the available phase space. We address 
this region specifically below. For the initial state wavefunction we use the ACCMM 
model. 

In this calculation we shall encounter two types of vertices: "penguin" and "tri- 
angle" . 

1. "Penguin" vertex (Fig. 2a). This vertex was calculated quite a while ago 0] 
and for the fost^-interactions turned out to be 

= -^VuV,,s^-l^^{Fl{x){e^, - k,m + 75) - m,F|(x)za,.r (l - 75)) (1) 

where x = mf/m^r for i = c,t. Fi{x),F2{x) are the Inami-Lim functions, Ff = 
Efc FgVbkVks, using unitarity of CKM matrix we get Ff = -5.24 and F^ = -0.19 for 
rrit = 176 GeV,mc = 1.3 GeV. Only the magnetic-moment type vertex contributes 
to the real photon emission in the b 57, so 

= ^ g^2 VbiVismbF^{x)ia^uq''{l - 75) (2) 

F2{x) ~ 0.65. We refer to k being the momentum of the gluon and q - momentum 
of the photon. In the gluon vertex, we shall drop the small contribution from , 
which, however, might be sizable for gluonic penguin-mediated hadronic processes 

i- 

2. "Triangle" vertex (Fig. 2b). The calculation of this vertex reduces to the cal- 
culation of the set of usual triangle diagrams with certain form-factors ( a complete 
parameterization of the diagrams of this type was first given in ||^) and is interesting 
in the sense that light (spectator quark) current couples to b —>■ s current through 
e-tensors. The vertex itself was recently calculated by D.Wyler and H.Simma in ||^ 
and takes the form 

= - ^ ^g^^ y Vfr^'^^s {ifi iius + if2 iius + i/s m^s)7'(1 + 75) (3) 

where 

/l tins ^rtivs{k ~\~ q ^^g) (^r^ius^ 

f2 ,us = '-^k'^q^k^M^, (4) 
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with 

A^5 = -1 + -^^{Qoizo) - Qo(^i)) ^(Q-(^o) - Q-i^i)) 

Zq — Zl Zq — Z\ 

A^e = 1 + ^^(Qo(^o) - Qo(^i)) + (Q-(^o) - Q-(^i)) 

Zq — Zl Zo — Zl 

Ai23 = Ais = -Ai26 (5) 

/■I 

= / — ln{l — ^^(l — u)) 
Jo u 

Qq{x) = du ln{l — xu{l — u)) 
Jo 

and Zq = s/m],zi = k'^/mf, s is a total momentum of internal quarks. 

These vertices will be used to calculate the rate and spectrum in the next section. 
We define the total decay rate as 



tot 



~l~ Tfic ~l~ ^mix (6) 



where Tg represents the leading contributions from b ^ s'y (diagrams of the class 
J), Tsc gives the corrections which involve the spectator quarks, and Tmix gives the 
mixing of the two. In powers of the strong coupling constant, Tmix is of order 
and Tsc is of order a^. 

In order to account for bound state effects, we will treat the decay of a i?-meson 
as a 2 — 3 scattering process, assuming that the scattering occurs from the bound 
state and model this state with the ACCMM- model @]. In calculating the decay 
rate (cross section 2 — 3) we average over the spins of the initial particles which, 
in principle, leads to the prediction for Ttot for the "mixture" of B and -B*-particles 
(one can, at the expense of a considerable increase in complexity of the calculation, 
separate B and B* by introducing polarization operators for projecting out different 
spin combinations). Thus the decay rate is defined as 

dT = - ^— \Tfi\^ (7) 

where factor of 4 in the denominator comes from the spin averaging and 

d^ = {2TiY5\pf-p,) , '^ll' , ^\t^ , (8) 
V 1 \yi -^^^ (2vr)32E, (27r)32£;4 (27r)32w ^ ^ 

We also treat spectator quark as massless to simplify the following calculations. The 
final phase space integration is performed numerically. 
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3 Spectator Contributions to the Rate and the 
Shape of the Photon Spectrum. 

3.1 Corrections of order 

The introduction of interactions with the hght quark obviously makes the kinematics 
of the process more involved - the final state is no longer a simple two-body final 
state which results in delta-function-like peak in the photon energy distribution. 
Instead, the peak is now extended down to lower photon energies making it possible 
for the photon to populate the energy interval from to ~ 2.5 GeV. Moreover, a 
difference between the decay spectra of neutral and charged mesons also arises. 

The question that naturally appears here is whether it is possible to have an 
interference of the spectator and non-spectator mechanisms. Since the 6 ^ 57 
transition is of order and the non-spectator matrix elements are first order in 
as this would provide corrections of order as- If there were not any momentum 
associated with the bound state, the spectator quark would remain precisely at rest 
in the 6 — > 57 transition and therefore would not interfere with the diagrams where 
momentum was transferred to the spectator quark. However if bound state effects 
are present the spectator quark will always carry some momentum (typically of order 
1 GeV or less) and the two classes of diagrams can interfere. Since this occurs only 
over one portion of the phase space, we can say that the interference effects are 
really of order aglGeV/mB- 

In the bound state, the initial momenta are described by a momentum distribu- 
tion (which is given by the wave function of quarks inside of the meson). The 
interference takes place with a probability proportional to a wave-function overlap. 
We now turn to the calculation of this overlap. 

Let Oi to be an operator responsible for the transitions via mechanism / and 
On - via mechanisms // — IV. Then, the mixing term arises as a result of 

\Mfi\^ = Y.{B'\Oi + On\Xs){Xs\Ot + 0+ |B) = 

= {B'\{Oj + Ou){Oi + Oij)\B) = {B'\\Oj\^ + \0n\' + 2Re OiO+\B) (9) 

where we have used completeness relation for |Xs)-states. We are interested in the 
last term. 

Let us write down the contributions from different classes of diagrams in turn. 
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^ ^ 2puq 2p'^q ^ > 

Here e is the polarization vector of the photon Qb, Qsp are the charges of b-quark and 
spectator respectively, and is given by (|l]). Note that fc^ in the denominator due 
to gluon propagator is to be cancelled by in fes^f-vertex, so Tn is a local operator. 
This operator (convoluted with the lowest order non-spectator graph, along with its 
charge conjugate part) gives a contribution to the shape of the photon spectrum. 
The result of convolution can be represented as 



2psq Pbq ^ 



Q 

vn^n r)i.n \ / 



KPb)lul9(Tap{"^ - l^>)Kp'b) + TTT^ Ps/ KPb)lt,lu(yap{"^ " 75)^(^6) 

^PuQ Puq 

{Su{qf,gae - Qfiaqe) - iPu(^f^atuqt ^5)} + C.c. 

with Sb = Psq+Ptq, Pb = Pbq-Psq and analogously for 5'„ and Pu, = u{p'jYu{Pu), 
L5 = u{p'^)'y^'y5u{pu) , = 4/3 is a color factor coming from the averaging over 
initial colors and summing over the final ones (taking into account the fact that 
initial and final states are color-singlets). Also 

_ I Ps P'u 

Pbq Puq Psq p'uQ 

Hpbht^luCTafiil - l5)b{p'i,) = {g^^vQaP + 
guag^l|3 - g,,ag,y(S - ie^^o/3)6(p'J(l - l^)h{pb) - (12) 

ib{p'b){g,^ucrap + guaCTfi/s - gtMaCTiyfs - iei,uae<ye(5){)~ - l5)b{Pb) 
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The calculation of the contribution from the diagrams of the class /// is more in- 
volved mainly because of the complicated structure of the vertex (Q) . The amplitude 
squared reads 

Th = Y^HfL^W., = -^h-F (13) 

with 

Const = -^m,\VtsVtb\^aQ''asClF:^{x) 

F = u{p'^)ru{Pu) (14) 
h, = 2b{p') ^ Uo.saauq'^l'il + iMp) 

i 

We will be working explicitly with components of this current in order to get iTjjjl"^. 
Using various 7-matrix identities 

i r 

lalu'ye = Qaule + Quela - Qaelu - i(^ave(3l(3l^ (15) 



we arrive at 

= 2zgVKp') H ft^as{9ue7a7s - QaBlvls + ieauBp'lpls]{l + l^)h{p) (16) 

i 

Note that the term symmetric under a ^ v disappears. Next, using explicit form 
of fifius, identities 

a'^^Ts = —e^'^-ax. (17) 

to get rid of e-couplings and discarding terms proportional to (on-shell photons), 
we write 

Ka = Hp'){p^s{{(11) + 2A^23(fcg)) - P(p.g) + q^{2{pj) - A^t2^{p,k) + ^t2,{Psq) + 

+2Ai23e^) - P(2A^23(Psg)) + 2P(2A^23^cr"^A;„g^ + A^23cr"^gaPs;3)(l^ 
+2z(p,g)a^^7,, - 2z(g/)a^>,, + 2i{psl)a^'qa - 2zAz23(PsA;)a^"g„ + 2zA«23((Psg)cr^"A;„ 



-^a^^gj - ^A^23(A:g)a'^>.„ - a"^gj^a^>,,,}(l + ^,)b{p) 



Using (|I4D we obtain the contributions from class ///. Finally, class IV gives 
the following 



Let us also mention that some range of momentum transferred to the spectator 
gives a contribution to the decay rate which is highly suppressed - in the ACCMM- 
model, for instance, decay rate looks like 

r = / d^vW)Y^^<^s{v) (20) 

withjf = p+fc and 0(j9) = {l/p^^)exp{—p'^/2p'^^). Note that this expression reduces 
to the conventional ACCMM-averaging if the momentum flow between heavy quark 
and spectator lines is discarded. It is clear that if the magnitude of 3-momentum 
transferred to the spectator line is large, the factor of wave function suppression is 
large as well which results in the exponential fall-off of the a^-spectrum as i?^ ^ 0. 
This is the reason why non-spectator corrections turn out to give much smaller 
impact on the line shape of the spectrum than those at order a^, in fact, 

Trm. ^ _ ^Q_2 (21) 



r. 



c 



3.2 Corrections of order a]. 

Let us now consider corrections of order a^. In fact they can simply be obtained by 
squaring the sum of Tj presented in the previous section. We shall present only few 
of them emphasizing their possible experimental signatures and highlighting some 
steps of derivations. 

Class I. Diagrams of class / give contributions to and have been calculated 
elsewhere [Q. This contribution to the branching ratio Br varies from 2.0 ■ 10~^ to 
4.2- 10~^. The major uncertainties at present come from the choice of the renormal- 
ization scale in the effective operators and the KM elements, and so this branching 
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ratio can be predicted more accurately in the future as the KM elements are more 
precisely measured. 

Class II. Diagrams of this class give a new contribution to Ttot which comes 
from the emission of the photon from the spectator quark. 

The matrix elements for diagrams (6) and (7) ( Fig. 1 ) can be rewritten as 

Tf^ = liH.K^'' (22) 

and are given by (|10]). Taking the hermitian conjugate and multiplying by ( |10|) gives 
|Tjjp for which we can write 

\Tf^? = H,,L^:pW-^ (23) 

where 

pol 

H,, = Tr{A ^ r^} (24) 

where as before we set m, = and is a color factor. To obtain we have used 
the fact that 5-meson is a superposition of the color-singlet states, so 

Cl = Itrf^^^^) = 1^ (25) 
° 3^2222^9 ^ ^ 

Evaluating the traces we obtain for if^i^: 

H^u = ^^a.slVbiVisl'^Ff '^k'^[psf,Pbu + PsuPbfi - {PsPb)gtiu + if^tiuaf^PsPb} (26) 

where the last term vanishes upon contraction with L'^p for which we can write down 
(taking into account contraction with Wa/s)'- 

Lz = Q%{2{puPi){{pM + p:yu - g'-'ip'uPu)) - p'm'p: + pup'u - 

g'^iPuPu)) + 2p^p^(pX) + 9'''WuP2){PuP2) - pUpup'u) + (27) 

pMpu + p^p'J)) - 2p'u'P2iPuP2) - 2p:pi^{pM) 



9 



for the diagrams (6) and (7). These diagrams might be interesting from the experi- 
mental point of view since they give rise to different contributions from charged and 
neutral 5-mesons (the amplitude is proportional to the charge of spectator Qsp), so 
one, in principle, can observe this charge asymmetry 

However the small size of our prediction means that it will be difficult to make this 
measurement. 

Class III. Calculation of the diagrams of this class involves the effective vertex 
(^ for the sum of two diagrams - 8 and 9. As before, the amplitude squared can be 
written in the form 

iT-i 1 2 a uafiab T P-^M[l<^0 

H^L = Trfr^'^ rf A} (29) 
Expanding T^-^ as T"-^ = E^ ft^.asls{l + 75) we get 

ik 

= 8 ^ fiastifkatu{pssPbt + PstPbs " {PsPb)gst + stnmP'^ PT} (30) 
ik 

Thus, the problem is the calculation of J2ik fias^ifkatu- This calculation, which re- 
quires the explicit form for the /-functions and involves expansions of the products 
of e-tensors, results in a quite lengthy final expression. However, restricting our 
attention only to the case of real photons (g^ = 0) and taking advantage of the 
symmetry of L^'^ under fi ^ u (implying further convolution with "light" tensor 
L^'^ ) simplify the result enormously since there are only four combinations of /'s 
that give non- vanishing contributions to the problem of interest. They are 

/l fj.saflusl3 = —{I'^ig^ivQalS — gavQ^ip) + ^/illSgau + J-uJ-agiM/S — ^al/Sgiiu — ^i^vgst) 
t f \ f f \ f f — /^^23^,2 I r,^^23A^6\ ,p aim n 



h MS"/3 usp = -q^qu^^^Y {{kq){qakp + qiiK - {kq)gap) - gag/3fc2}(31) 
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iisafl us[3 



1 iMsaj3 vsP 



{lk)g^p) - qpq^{kJa - {lk)gav) + Qapik^q^ + q^k^){lq) - {lq){kaqu9,ip + k^q^gau)^ 

Combining all of the terms together and contracting iJ^" to L^^ one obtains the 
contribution of the class III diagrams to the Ttot which is to be computed numerically. 
The resulting analytical expression for the amplitude is too long to be presented here. 

Class IV. Finally for the class of diagrams responsible for the 6(s)-spectator 
hard gluon exchange, using the same technique as above, we write for L and H 
tensors 



jflV 



Here the "heavy" current is 

1 



H^at = Tr{^s n;: A n^} 



^1 



(32) 



(33) 



and is given by (|^). Evaluating the traces and making use of the facts that L 
is symmetric under fi u and = 0, we obtain 



ab 



TjaOL 



2^ 



32 



k'M 



rri] 



2\2 



P\ - ^1) 



{pU^ + pW - {Pbq)9n - 2{pqMpl + p^ - {PiPb)9n) + (34) 

^Mq){pMq'' + pW - {Psq)gn - 2(p2g)(p>^ + - {PsP2)gn)_ 



with pi = ps + q and p2 = Pb — q, here Qint is a charge of the internal quarks. 
Contracting -ff^^^fe with L'^'j^ gives the expression for the diagrams of class IV. In the 
final calculation of the magnitude of the spectator effects, the dominant contribution 
comes from the class IV diagrams 0. 

In evaluating the matrix elements, we have taken the 6-quark spinors to be non- 
relativistic. The interferences between the different classes of diagrams involving the 
spectator quark have also been taken into account, although we have not displayed 
the (lengthy) formulas explicitly above. Some of the complicated sums over many 
Dirac matrices were evaluated numerically with the help of Mathematical^' . 
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Substituting \Tfi\^ and ( ^HD into (|^, and performing numerical 5-dimensional 
Monte-Carlo integration over (j)s,(l>q,9q,9s,p we obtain photon energy distribution 
presented in Fig. 3. Integrating over x = 2uj/Mi, gives for the decay rate: 

6Ttot ~ 5.3 ■ 10"^^ GeV (35) 
55r~1.5-10-^ (36) 

Given the recent measured by CLEO value for Br{B Xgj) = 2.32 ■ 10^^ 
we conclude that spectator corrections give an overall effect of about 5% on the 
inclusive decay rate T{B — >• Xsj). 



4 Summary 

The diagrams which involve the spectator quark are in general non-local and involve 
momentum transfers to the spectator which are comparable to the energy released 
in the decay. This is a quite different configuration from the 6 — > 57 transition 
where the spectator quark is nearly at rest. Despite the inclusion of the bound 
state momentum of the spectator, the interference of the two classes of diagrams 
turns out to be very small. We have evaluated the remaining effects of the spectator 
diagrams, which then amount to effects at order a^. 

Our results indicate that the diagrams involving the spectator quark do not play 
a large role in 5 — »• ^^7, which implies that the 6 — > 57 reaction is the dominant 
contribution to the inclusive decay. 

Acknowledgments. One of us (A. P.) would like to thank S. Nikolaev for useful 
conversations regarding Monte-Carlo integration. This work has been supported in 
part by the US National Science Foundation. 

Appendix. Calculation of the Phase Space. Since we are calculating 
Q CD- correct ions associated with bound-state effects it is convenient to work in the 
rest frame of B-meson. Let us adopt ACCMM-model for use in this particular 
frame. In this coordinate system (the center of mass frame for the system "heavy 
quark-spectator" ) : 

^u + Ps = = Pu+Pb (37) 
Eb + Eu = Mb = E'^ + Es + iu (38) 

We treat the spectator massless particle, so 

Eu=Pu= P, Eh = ymf+p^ (39) 
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As usual in the ACCMM model the effective b quark mass is momentum dependent 
and is given by ml = Mb{Mb — 2p). Also 



Es = VpI = yiP'u + q)' = ^jE'i + + 2E'^ujcos{e'J (40) 

Let us note that in this frame the decay products m', s, q lay in the same plane, 
thus the angle parameters are ^^^-the angle between the directions of 3-momenta 
of s-quark and photon and the relative angle between 3-momenta of the initial 
( "spectator" ) w-quark and normal to the decay plane. Let us recall that 

5\Pf - Pi) = SiEj - Ei)5\pf - Pi) (41) 
Therefore, (||) can be integrated over to get 

^ ^^5{Ms -E:-E,~uj) (42) 
uj E^Es 

with = {Eh + Eu)"^. Next, evaluating the p'^-integral (keeping in mind that Es 
is given by (|40|)) yields 

d^qdQ'^E'udE' , 

d{MB - E^- Es - UJ) (43) 

u Es 

Making use of a familiar ^-function relation 

^(/(^))= E n77^^(^-^o) (44) 
to perform integration over E'^ we arrive at the equation for the roots of E'^: 



Mb-E'^- ^E'l + uj^ + 2E'^u:cose'^^ - u; = (45) 
which can be solved to get 



Mb{Mb - 2u) 



2{MB-uo{l-cose')) 



(46) 



This gives for the (5-function: 



Hence, the phase volume is 



" "2(MB-^(l-coseg)2 ^ 



^yiv±B — — ^^■^^uq)^'^ 

with dVLi = dcosO, 
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Figure 1: Set of the relevant diagrams. 
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Figure 2: Effective vertices. 
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Figure 3: Spectrum for B Xs'j. 
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